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Spectral positivity is known to be violated by some forms of lattice gauge fixing. The most notable 
example is lattice Landau gauge, where the effective gluon mass is observed to rise rather than fall 
with increasing distance. We trace this violation to the use of quenched auxiliary fields in the lattice 
gauge fixing process, and show that violation of spectral positivity is a general feature of quenching. 
We illustrate this with a simple quenched mass-mixing model in continuum field theory, and with 
a quenched form of the Ising model. For lattice gauge fixing associated with Abelian projection 
and lattice Landau gauge, we show that spectral positivity is violated by processes similar to those 
found in quenched QCD. For covariant gauges parametrized by a gauge-fixing parameter a, the 
517(2) gluon propagator is well described by a simple quenched mass-mixing formula. The gluon 
mass parameter appears to be independent of a for sufficiently large a. 

m . 

O ; 

Although many observables can be determined in lattice gauge theories without gauge fixing, there are several 
reasons why gauge fixing is desirable in lattice simulations. Gauge fixing is necessary to make the connection between 
continuum and lattice gauge fields. Continuum theories of the origin of confinement often make predictions about 
the gauge field propagator. Gauge fixing has also been a key technique in lattice studies of confinement as well.Q 
Important properties of the quark-gluon plasma phase of QCD, such as screening masses, are contained in the finite- 
temperature gluon propagator. 

Techniques for lattice gauge fixing have been known for some time. It has been clear from the beginning that 
■ non-Abelian lattice gauge field propagators show a violation of spectral positivity. This is readily seen from the 
' effective mass: for a normal operator which connects only states of positive norm to the vacuum, the effective mass 
monotonically decreases with distance to the lightest mass state coupling to the operator. Covariant gauge gluon 
propagators have an effective mass increasing with distance. In one sense, this is not surprising. We know from 
perturbation theory that covariant gauges contain states of negative norm. However, that knowledge has neither 
explained the form of the lattice gluon propagator nor aided in the interpretation of the mass parameters measured 
(f) . from it. In fact, no similar violation of spectral positivity is observed in the U(l) case Q, which has negative-norm 
<0 states in covariant gauges. 

In lattice simulations, gauge fixing has typically involved choosing a particular configuration on each gauge orbit. A 
brief review of this approach is given in Ref. |J] . In the continuum, on the other hand, gauge fixing usually includes a 
parameter that causes the functional integral to peak around a particular configuration on the gauge orbit. As shown 
below, the extension of this idea to lattice gauge theories makes clear that lattice gauge fixing is a form of quenching, 
. with the gauge transformations acting as quenched fields. As has been demonstrated in quenched QCD, quenching 
can violate spectral positivity, with significant effects on many observables. 0, 0] 
• »"H . We begin with a review of lattice gauge fixing, including the generalization of lattice Landau gauge to covariant 
gauges with a gauge parameter. This generalization will be directly interpreted as a quenched Higgs theory. We 

then explore the origin of violations of spectral positivity in some simple lattice and continuum models of quenching. 
Simulation results for the effective mass of an SU(2) lattice gauge field will show behavior very similar to these models 
as the gauge fixing parameter is varied. We will argue that spectral positivity violations in both lattice covariant 
gauges and in studies of Abelian projection originate in the quenching process. 

The standard approach to lattice gauge fixing is a two step process. |^ An ensemble of lattice gauge field configu- 
rations is generated using standard Monte Carlo methods, corresponding to a functional integral 
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where Sjj is a gauge-invariant action for the gauge fields, e.g., the Wilson action. The gauge action Sjj is invariant 
under gauge transformations of the form (x) — ► g (x) (x) g + (x + /i) . 

In order to measure gauge-variant observables, each field configuration in the [/-ensemble may be placed in a 
particular gauge, i.e., a gauge transformation is applied to each configuration in the [/-ensemble which moves the 
configuration along the gauge orbit to a gauge-equivalent configuration satisfying a lattice gauge fixing condition. The 
simplest gauge choice is defined by maximizing J2 X M Tr [[/ M (x) + J7+ (x)] for each configuration over the class of all 
gauge transformations. Any local extremum of this functional satisfies a lattice form of the Landau gauge condition: 

Y,lAA* + ri-^(x)]=0 (2) 
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where (x) is a lattice approximation to the continuum gauge field, given by 



U,(x)-U+(x) 1 
M (X) = 2z 
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Other gauge-fixing conditions may also be used ficl ]. and lattice improvement techniques can be applied to the 
definition of to reduce discretization errors as well. The global maximization needed is often implemented as a 
local iterative maximization. The issue of Gribov copies arises in lattice gauge fixing because such a local algorithm 
tends to find local maxima of the gauge-fixing functional. There are variations on the basic algorithm that ensure a 
unique choice from among local maxima. |lfl| 

For analytical purposes, it is necessary to generalize this procedure ,9|, so that a given single configuration of gauge 
fields will be associated with an ensemble of configurations of g-fields. We will generate this ensemble using 

S 9 f [U, g] = ^ Tr [9 ( x ) U » (*) 9 + & + M) + .9 (* + /x) U+ (x) g+ (x)] (4) 
I 

as a weight function to select an ensemble of g-fields. The sum over I is a sum over all links of the lattice. The 
normal gauge- fixing procedure is formally regained in the limit a — > oo. Computationally, this can be implemented 
as a Monte Carlo simulation inside a Monte Carlo simulation. 

Note that the ^-fields must be thought of as quenched variables, since they do not affect the [/-ensemble. The 
expectation value of an observable O, gauge-invariant or not, is given by 

(O) = 4~ I \- dU ^ eSu[U] ^TF^ i \ d 9\ e S ° f[u ' 3] O , (5) 
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where 



[[/] = J [dg]e s ^ u ' g \ (6) 



Formally, the field g is a quenched scalar field with two independent symmetry groups, G g i b a i ® Gi OC ai, so that it 
appears to be in the adjoint representation of the gauge group, but the left and right symmetries are distinct. The 
generating functional Z g f[U] is in some ways a lattice analog of the inverse of the Fadeev-Popov determinant . 1 1 lj 
However, there are important differences. Note immediately that Z g f[U] depends on the gauge-fixing parameter a. 
More fundamentally, the lattice formalism resolves the Gribov ambiguity. By construction, gauge-invariant observables 
are evaluated by integrating over all configurations. Gauge- variant quantities receive contributions from Gribov copies, 
always with positive weight. Thus the connection between this formalism for lattice gauge fixing and gauge fixing in 
the continuum is not simple. Furthermore, alternative lattice gauge fixing procedures have been proposed,along with 
new gauge choices specific to the lattice. A comprehensive review of gauge fixing technology is available. [lOj 

We begin our analysis of spectral positivity violation with the simplest model of quenching possible: two free, real 
scalar fields with a non-diagonal mass matrix. The Lagrangrian is 
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We treat the quenched approximation of this model in a manner completely parallel to our discussion of lattice gauge 
fixing above. We divide the action into three parts S = S\ + S2 + Si 2, where Si and S2 are functionals only of 4>\ 
and 02, respectively, and Si 2 contains the mixing term. We quench the field <f>2- Although there are no loops in 
this simple theory, quenching implies that 02 cannot appear as an internal line in the complete propagators. The 
generating functional in the quenched approximation, including sources J± and J2 is 
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where we have introduced a kind of ghost variable 02; space-time variables are implicit. 
From the generating functional we can obtain the (0i0i) and (0202) propagators. 
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In momentum space, the 



f>i0i) propagator is X/(p + mf), since 0i is unaffected by 02- On the other hand, the (0202) propagator is 



1 



1 



1 



(9) 



3 



+ ■ ■ 

FIG. 1: Exact propagator for the quenched <j>2 field in the simple mass mixing model. 

An alternative diagrammatic procedure is to sum Dyson's series, as shown in Fig.Q noting that the ((^2^2) propagator 
is truncated at two terms. The propagator has a structure similar to the 7/ propagator in quenched QCD the 
77' has a double pole form in quenched QCD when singlet self-energy graphs are approximated by a constant. The 
(02 ^2) propagator also may be written as 

L s \ 1 , s 1 , s 1 (10) 



\ (m 2 2 - mlf J P 2 + ml (m 2 - m 2 ) 2 P 2 + rn 2 mj - m 2 2 ( p 2 + m 2f 

This propagator always violates spectral positivity because of the double pole term, 1 / (p 2 + tttJ) 2 , which has a 
coefficient whose sign depends on m\ — rn 2 ,. Another possible violation of spectral positivity occurs for sufficiently 
strong mixing: if /i 4 > (ml — mf) , there is a simple pole at p 2 = —ml with negative residue. 

The form of the (</>2<^2) propagator in coordinate space is very interesting, and forms the basis for our study of other 
quenched theories. In any number of dimensions, we can consider propagators using wall sources, i.e., of co-dimension 
1. This has the effect of setting the momentum equal to zero in all the directions of the wall. For wall sources, we 
have the propagator 
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The factor 7772 \x\ e m2 l x l shows an initial rise rather than a decay with increasing violating spectral positivity. 
We define an effective mass associated with the <p 2 field as 

m eff (x) = lim /in ( f [X } ) = -±ln(G(x)). (12) 
a^o a \G(x + a) J ax 

One can easily check explicitly that m e ff(x) — > min(mi, 777,2) as a; — > 00. For any field theory which obeys spectral 
positivity, m e ff(x) monotonically decreases to its limiting value. Theories violating spectral positivity may display a 
complicated behavior for m e ff(x) before the eventual onset of asymptotic behavior. 

We have identified three different possible behaviors for m e ff(x) in this simple quenched model. If the mixing 
parameter \x is sufficiently small and m\ < m 2 , m e ff(x) monotonically decreases to its value at infinity, as in a 
normal field theory which obeys spectral positivity, as shown in Fig. EJa). As fi is increased relative to mi and 777.2, 
m e ff(x) may develop a minimum, as displayed in Fig. EJb). On the other hand, if 777,2 < mj, the behavior seen in 
Fig- 121a) is not possible, and only the behaviors seen in Fig. Etc) and Fig.|2td) are possible. In Fig. Etd) , the minimum 
has moved to x = 0. For sufficiently small /i, these effects are difficult to observe, and m e ff(x) is essentially equal to 
77Z2 for all x. Regardless of the relative size of mi and m 2 , an observable violation of spectral positivity associated 
with m e ff(x) not monotonically decreasing indicates a significant mixing parameter \x. 

Similar behavior can be observed in a very simple lattice model based on the Ising model, where real-space arguments 
can be used to find an approximate propagator. We consider two coupled one-dimensional Ising models, with spins 
fii, <Ji e {— 1, +1} and respective nearest-neighbor couplings J and K. The a spins are coupled to the /j spins via 
an interaction of the form ^ i Lai/ii, and the cr's are quenched. This simple model is a form of spin glass, with the 
averaging over the ensemble of /i spins representing the "quenching" process. 

The a propagator is given by 



-7TT~\ Yl 170(7 n exp 
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where is the partition function for /i and Z a [(i] is the partition function for a in the presence of a particular fi 
background. The parameter L is a mixing parameter. We can approximately evaluate the a propagator for J, K, and 




FIG. 2: The effective mass associated with <f>2 as a function of x. The parameters used are: (a) mi = 0.1, 7712 = 0.4, /i = 0.2; 
(b) mi = 0.1, m 2 = 0.2, fj, = 0.2; (c) mi = 0.5, m 2 = 0.4, fi = 0.4; (d) mi = 0.5, m 2 = 0.1, fx = 0.4. 



L sufficiently small by considering the direct contribution (tanh if)™ combined with mixing of crwith [i. This indirect 
term can be written as (compare Fig. ^) 

n n—p 
p— 1 m— 

After performing the summations, the propagator is given approximately as 

(cro<T n ) w (tanhX)™ + (tanhL) 2 (tanh if) 

where a; = tanh J/ tanh if. The n(tanhi4T)™ factor signals a violation of spectral positivity, just as the 
m \x\ exp (— m \x\) term did in the mixing model. Of course, the arguments which led to Eq. © and Eq. 115|) 
are essentially the same, but carried out in momentum space and real space, respectively. For small J, K, and L, 
Eq. i|15|) fits lattice simulations of the (<J$<J n ) propagator well. 

In Fig. 13 we show the effective mass determined from the {uo<y n ) and (^oM«) propagators for the parameter set 
,/ = 0.7, K = 0.9, and L = 0.3 for a one-dimensional lattice of size 26. The propagators were obtained from 40000 
heat bath sweeps of the fj, variables; after each such sweep, 100 heat bath sweeps of the a variables were carried out. 
The parameters K and L were chosen empirically so as to display a clear violation of spectral positivity. The fj, mass 
fits very well with the analytical solution m = — lntanh(J) for the d — 1 Ising model out to a distance of 8. Note 
that the a reaches its asymptotic value of — ln(tanh_ftT) from below, and only at n ~ 8. The similarity to the simple 
field theoretic model of quenching is clear. 
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FIG. 3: Effective masses for a and n in the 1 — d Ising model. 



We will now show that the SU{2) lattice gluon propagator regarded as a function of a shows behavior similar 
to that of the other, simpler quenched models studied above. Simulations of this type of lattice field theory, with 
stochastic quenched gauge fixing fields, were first performed by Henty et al. [l2j, who studied the case of SU(i) as a 
function of a at /3 = 5.7 on 8 4 lattices. They found evidence for a first-order phase transition as a was varied, but did 
not determine the full phase diagram in the a-(3 plane. They also found that the gluon propagator was dependent on 
a, a result which could be anticipated from the strong-coupling expansion. 

Let us consider for the moment the unquenched version of the gauge fixing model. This is a model with scalar 
fields in the fundamental representation of the gauge group in addition to the gauge fields. The scalar fields explicitly 
break the Z{N) global symmetry associated with confinement in the pure gauge case, and external color charges 
are screened. As first shown by Fradkin and Shenker this leads to a connection between the strong-coupling, 
confining phase and the Higgs phase, so the two phases are not actually distinct. We have verified that this phase 
structure is preserved in the quenched form of the model. For (3 sufficiently large, there is a line of first-order phase 
transitions in the (3-a plane. It is very reasonable that such a line exists in the quenched model, since it can be 
thought of as the continuation of the critical point of a pure spin model at (3 = oo. However, this line terminates 
at a critical end point; for sufficiently small a, the nominal confining phase (/3 small) and Higgs pha se {(3 large) are 
directly connected. This observation forms the starting point for a detailed analysis of the model. [l4j 

We have performed simulations of SU(2) gauge theory at j3 = 2.6 with a ranging from 1.0 to 3.0 on a 12 3 x 16 
lattice. At this value of /?, there is a first-order phase transition at a 0.83. We have fit the data using a simple 
generalization of the quenched mixing model. The coordinate space propagator has the form 



(A + Bm 2 \x\)e~ 
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This form for the propagator follows from the replacement of the mixing parameter [i in Eq. @ by the more general 
form /j 2 (p 2 + m 2 ) . 

In Fig. 21 we show the effective mass as a function of x for a — 1, 2, and 3. The solid lines are obtained from fits 
to the propagator using Eq. i|16|) . The similarity to the other quenched models is quite clear. For small a there is an 
initial decrease and then rise of the effective mass, much like Fig. |2Jc); as a increases, this minimum vanishes. We 
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FIG. 4: Effective masses for three values of the gauge fixing parameter, a — 1.0, 2.0 and 3.0. 



plot in Fig. [3] the best-fit values of mi and m-i as a function of a. Note that the behavior of mi is consistent with it 
being constant in this region, while m± appears to be decreasing to a constant limit as a increases. 

Our results are roughly consistent with the work of Leinweber et al, who performed high-precision studies of the 
SU(3) gluon propagator at a = oo.^j| Among a large variety of possible functional forms for the gluon propagator, 
they found that their data was best fit by the functional form 



G(k) = Z 
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with L{k 2 , M) an infrared-regulated version of the asymptotic behavior of the renormalized gluon propagator in the 
continuum. Their best fit was achieved with the parameters 6 = 2.2to\l±° Q 2 3 , M = (1020 ± 100 ± 25) MeV, and 
A = 9.8^0 g- Many other functional forms were ruled out. 

Our results suggest that the lighter mass parameter m-i is independent of a, at least for large a (in the Higgs 
phase). If m 2 is indeed independent of gauge choice, as least within the class of covariant gauges considered, it seems 
natural to identify it as the gluon mass. As a consequence of the quenched character of lattice gauge fixing, this state 
partially mixes with another, heavier state, with a mass on the order of the scalar or vector glueball. 16] 

Note that the value of the lightest mass in the propagator may be difficult to extract from the effective mass. While 
it is true that the effective mass tends asymptotically to the lightest mass, the approach to the limit can be much 
slower than in a conventional field theory obeying spectral positivity. For example, at a = 3.0, m e jf at x — 7 is 
substantially lighter than m%. Having a theoretical basis for the form of the propagator is crucial in estimating the 
mass. 

Another application of lattice gauge fixing is Abelian projection, a method for investigating the confining properties 
of gauge theories. In lattice gauge theories, Abelian projection is implemented as an algorithm for extracting an 
ensemble of Abelian gauge field configurations from an ensemble of non- Abelian configurations. A notable success 
of lattice studies of Abelian projection [tH Hs| has been the correlation of the string tension of the projected theory 
with the string tension of the underlying non- Abelian theory. 




FIG. 5: Values of the two mass parameters in eq. 1161 as a function of a. The light mass, m,2 is approximately constant at the 
value of 0.283(7), while the heavy mass initially decreases with increasing a, reaching a constant value of 1.151(8). 



These lattice studies of Abelian projection typically use gauge fixing in an integral way. Taking for clarity the case 
of SU(2), the gauge fixing functional is 

S 9f = E f Tr [9 (z) U„ (x) g+ (x + /i) a 3 g (x + n) U+ (x) g+ (x) a 3 ] , (18) 

33, fl 

which is conventionally maximized over the gauge orbit, corresponding to the limit a — > oo in the formalism used 
here. An initial study of the phase structure in the SU(2) case finds evidence for a first order phase transition as a is 
varied at (i = 2.4. [ig The aim of this procedure is to transform an SU(2) configuration into a gauge-equivalent one 
which lies mostly in a given U(l) subgroup. After this gauge-fixing, the actual projection to U(l) is performed. 

In the case where no gauge fixing is done (a = 0), and only projection occurs, Faber et al. and Ogilvie HU 
have proved that the asy mptotic string tension measured in the projected and underlying theories are the same. 
Furthermore, Ogilvie [2l| has proven that this result should continue to hold for small a, under the assumption that 
the gauge fixing does not violate spectral positivity. However, the fact that the string tension evaluated using various 
forms of Abelian projection with gauge fixing is consistently slightly different from the actual non-Abelian string 
tension [52, suggests that a violation of spectral positivity may indeed be occurring. 

We identify the origin of this violation as the presence of a quenched scalar field. The case of SU{2) is particularly 
clear. Note that the combination g + (x) a^g (x) occurring in S g f can be written as a Hermitian scalar field (f>(x), 
where <f> transforms as the adjoint representation of the gauge group. The field <f) is traceless, Tr{4>) — 0, and satisfies 
Tr(</) 2 ) = 2. The gauge fixing action is thus equivalent to an adjoint scalar action of the form 

Sgf = E \ Tr y> W u » (*) * ( x + K (*)] • ( 19 ) 



As we have seen, such quenched fields naturally lead to violations of spectral positivity. Suppose we wish to measure 
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FIG. 6: An example of a problematic diagram when calculating the expectaion value of Eq. 1201 . The solid line is the Wilson 
Loop, the dotted lines are <j> propagators, and the wavy lines are gluons. 



a U(l) projected Wilson loop. This may be obtained from the expectation value of 

Tr J] 5(1 + cx z ) gj U jg j +x = Tr J[ 1(1 + ^)U, , (20) 
j 

where the product is ordered along a closed path labeled by the index j. The U(l) projected loop is represented in 
the underlying quenched Higgs theory as a sum of Wilson loops with all possible insertions of <j> at lattice sites on the 
path. When four or more <f> fields are inserted, problematic subdiagrams appear of the type shown in Fig. H3 Such 
terms lead to a violation of spectral positivity: there are no internal <f> loops in the quenched approximation, and an 
infinite set of diagrams occurring in the full, unquenched theory is omitted. This exactly parallels quenched QCD. 
This work was partially supported by the U.S. Department of Energy under grant number DE-FG02-91ER40628. 
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